Hierarchical organization of brain functional network during visual task 
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In this paper, the brain functional networks derived from high-resolution synchronous EEG time 
series during visual task are generated by calculating the phase synchronization among the time se- 
ries. The hierarchical modular organizations of these networks are systematically investigated by the 
fast Girvan-Newman algorithm. At the same time, the spatially adjacent electrodes (corresponding 
to EEG channels) are clustered into functional groups based on anatomical parcellation of brain 
cortex, and this clustering information are compared to that of the functional network. The results 
show that the modular architectures of brain functional network are in coincidence with that from 
the anatomical structures over different levels of hierarchy, which suggests that population of neu- 
rons performing the same function excite and inhibit in identical rhythms. The structure-function 
relationship further reveals that the correlations among EEG time series in the same functional 
group are much stronger than those in different ones and that the hierarchical organization of brain 
functional network may be a consequence of functional segmentation of brain cortex. 

PACS numbers: 87.19.1j, 87.19. le, 87.19.lt, 89.75.Fb 



I. INTRODUCTION 

Human brain, which is consisting of ten thousand mil- 
lion neurons and even more synapses, is perhaps the most 
complex system ever known. Benefiting from the devel- 
opment of brain anatomy since the nineteenth century, 
we now know that the neuronal elements of the brain con- 
stitute an extremely complex structural network, which 
supports a wide variety of cognitive functions and neu- 
ral activities [l|, H|. Recently, it has gained great in- 
terests of scientists to investigate functional connectiv- 
ity of brain based on complex network theory in which 
the brain can naturally be abstracted as a functional 
network. A brain functional network can be extracted 
based on functional MM (fMRI), electroencephalography 
(EEG), magnetoencephalography (MEG), or multielec- 
trode array (MEA) data, which record electric, magnetic, 
or other signals representing cognitive cortical activities 
of brain Q . Vertices of brain functional network derived 
from fMRI data describe anatomically localized regions 
of interest (ROIs) or voxels of fMRI image, whereas the 
vertices of those derived from EEG, MEG, or MEA data 
mimic the surface electrodes or sensors. The functional 
connectivity (or edge) between pairs of vertices is usually 
estimated using correlation between records of different 
vertices. 

It has been widely observed that the brain functional 
networks demonstrate properties such as small- worldness 
[3| and power- law degree distribution Q, which distin- 
guish brain functional network from regular and ran- 
dom networks [fjl-tTlj]. However, the small- worldness and 
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power-law degree distribution only represent the global 
properties of brain functional networks. To understand 
brain functional networks imposed by structural and 
functional constraints more comprehensively, other as- 
pects, which are expected to reflect both local and global 
organization of brain functional network (e.g. hierarchi- 
cal modular organization), need to be investigated. 

Hierarchical organization, also called community struc- 
ture and modular architecture, describes that some nodes 
in a network arc densely connected as groups, and these 
groups are only sparsely connected among themselves, 
which is a common phenomenon in diverse networks such 
as World Wide Web, scientist collaborating networks, ge- 
netic networks, protein-protein interaction networks, and 
financial networks [12h15| | . A large number of algorithms 
are developed to detect the hierarchical organizations of 
real networks [HI, El-Ill ( also 

see review work (2lj|). In- 
terestingly, several recent works have also found the hier- 
archical organization of brain functional network derived 
from resting-state fMRI data [HHH and epileptic MEG 
signals [HI, respectively. 

Since the large-scale neuronal networks of brain emerge 
from synchronized delta oscillations, the cortical EEG 
time series are able to distin guis h the ongoing from the 
evoked activities of brain (25j, 26]. In our work, the brain 
functional networks are derived from high-resolution syn- 
chronous EEG time series, which consist of 238 channels 
and are recorded during a cognitive task involved in vi- 
sual, judgment and motor functions of brain. In partic- 
ular, vertices correspond to surface electrodes (i.e. chan- 
nels), and edges are determined by correlation (degree 
of phase synchronization) of EEG time series from pairs 
of channels. In fact, the electrodes can also be clustered 
into functional groups by their spatial positions on the 
scalp and specifically priori knowledge of anatomical par- 
cellation of brain cortex. For example, electrodes near vi- 
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sual cortex area mainly represent visual function, and are 
therefore clustered together. The anatomical parcellation 
of brain cortex is performed according to Brodmann seg- 
mentation scheme [27j . The fast Girvan-Newman (GN) 
algorithm is first applied to analyze the hierarchical orga- 
nizations of brain functional networks. Strong evidence 
for the existence of modular architecture is found. Com- 
paring the clustering results from the brain functional 
network and the anatomical segmentation, we find sig- 
nificant coincidence of the modular architecture in the 
network derived from EEG data and anatomical organi- 
zation of the cortex over different levels of hierarchy. This 
result suggests that vertices are more tightly coupled in 
same functional cortex region than those belonging to 
different ones and that the patterns of neural activities 
of brain cortex are to a large extent determined by the 
anatomical modular architectures of the brain. 



II. MATERIALS AND METHODS 

A. Data acquisition 

The high-resolution EEG time series were syn- 
chronously recorded during visual task by using a large 
number of scalp electrodes (238-channels). Therefore, 
they had a high spatial and temporal resolution, which 
provide very helpful and detailed information of electrical 
activity of cortical surface. Specifically, the single-subject 
data set was recorded by A. Delorme et al. in the Swartz 
Center, UCSD, with the sampling rate 256 Hz using a 
Biosemi Active Two system [28| . Experiments were per- 
formed as follows: Filled white disks appeared briefly in- 
side one of five empty squares, and one of the five outlines 
was colored green to mark the square as a visual target, 
then the candidate made a motor response by pressing a 
mouse button with their right hand as quickly as possible 
whenever the filled white disk appeared at the attended 
location. It is noted that these locations were counterbal- 
anced across blocks in pseudorandom order. The more 
detailed description of experiment can be found in sev- 
eral previous works [29|,[3(|. In this data set, each visual 
target was represented by a synchronous EEG recording 
with five sections, and about three thousands trials of 
experiment were performed. In addition, there are 235 
channels of EEG used in our work by dropping 3 channels 
of EOG. 



B. Network construction based on Phase 
synchronization 

The concept of phase synchronization is introduced 
to study synchronization of coupled oscillators and has 
gained particular interests to investigate coupling among 
nonlinear complex systems (3fl433| . The phase of a 
real- value time series X(t) is defined using Gabor's ana- 



lytic signal approach [34| : 

V(t)=X(t)+iX h (t)=Ae^ t \ (I) 

where the imaginary part Xh(t) is the Hilbert transfor- 
mation of X(t). Hence, the degree of phase synchroniza- 
tion between time series Xi(t) and Xj(t) is evaluated by 
a bivariate phase coupling index [35j |: 

^HI~X> i( * (t) -* (t)) ll, ( 2 ) 

t=i 

which is in the interval [0,1]. If the phases of two time 
series are completely synchronized, the phase coupling 
index will be maximum. 

The phase synchronization method is able to reduce 
the nonstationary effect, compared with the calculation 
of the correlation of two EEG time series directly. The 
functional connectivity (edge) between pairs of electrodes 
(vertices) has been estimated using the measure of phase 
synchronization, and then it is thresholded to generate 
functional networks. Herein each vertex is connected to 
its N nearest neighbors (i.e. those channels that are 
most phase-synchronized with itself) [36j. Notice that 
the resulting brain functional networks are not necessar- 
ily symmetric, i.e. if vertex i is a neighbor of vertex j, 
vertex j may not necessarily be a neighbor of vertex i, 
and vice versa. Though the edges are intrinsically di- 
rected, we analyze the networks in an undirected way for 
convenience. The mean degree of network is generally 
determined by N. However, since the directed edges are 
identified as undirected ones, the actual mean degree is 
a little larger than N. In this way, the connectivity of 
brain functional network is naturally guaranteed with- 
out dense connections, which will otherwise disturb the 
hierarchy detection. 

C. Clustering EEG channels by electrode position 
and anatomical parcellation of brain cortex 

EEG time series recorded through each channel mainly 
represent the electrical activities of neurons near the cor- 
responding electrode. Thus, the EEG channels can be 
directly clustered into functional groups according to the 
spatial position of electrodes and anatomical parcella- 
tion of brain cortex based on Brodmann area. We used 
the Brodmann template image distributed with MRIcro , 
which is restricted to the standard MNI space [37[ • Res- 
olution of image is 181 x 217 x 181 and size of voxels is 
lmm x x 1mm x 1mm. Each hemisphere is partitioned 
into 41 areas according to cytoarchitecture of neurons, 
together with the same labels suggesting same cognitive 
functions in two hemispheres. The locations of electrodes 
are registered to MNI space using SPM8 toolbox (open- 
source software) (38|. By checking Brodmann areas in 
which the electrodes belong to, the electrodes that are 
spatially nearby are labeled by the same Brodmann ar- 
eas index. At last, these 235 channels are clustered into 
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TABLE I: Number of channels restricted to Brodmann area. 



Brodmann 


EEG channel 


Brodmann 


EEG channel 
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number 
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number 
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FIG. 1: (Color online) Average cluster coefficient and charac- 
ter path length as a function of N 



TABLE II: Number of channels restricted to functional sub- 
structure 



Functional 
substructure 


Brodmann 

area 


EEG channel 
number 


Broca's(B) 


44, 45 


17 


Audition(A) 


22 


5 


Cognition(C) 


9, 10, 46, 47 


31 


Emotion(E) 


38 


5 


Vision(V) 


17, 18, 19 


50 


Vision-parietal(Vp) 


7, 39 


18 


Vision-temporal(Vt) 


20, 21, 37 


51 


Motor(M) 


4, 6, 8 


31 


Sensory (S) 


1, 2, 5, 40 


19 




Number of neighbors 



25 functional groups corresponding to 25 Brodmann ar- 
eas (e.g. functional group 47 corresponds Brodmann area 
47). Other Brodmann ares are absent due to deeply in- 
side position of the brain cortex. The size of functional 
groups varies from 1 to 29, as shown in Tab|I] Moreover, 
Brodmann areas can be roughly scaled into 9 major sub- 
structures according to specific physiological functions, 
by which the EEG channels can be further clustered at 
a higher hierarchy, as shown in Tab[II] Hence, the EEG 
channels are organized into two levels of hierarchy ac- 
cording to the above mentioned segmentation. 

The hierarchical organization of EEG channels have 
two aspects, modular architecture of brain functional 
networks detected by fast GN algorithm and functional 
groups according to anatomical parcellation of brain cor- 
tex based on Brodmann area. Thus, the association of 
these two kinds of hierarchical organizations, e.g. the 
overlap between the communities of the channel sets and 
the functional groups, is a key question which is expected 
to shed lights onto the relation between brain function 
and structure. 



FIG. 2: (Color online) Maximum modularity and correspond- 
ing community number as a function of N 



III. EMPIRICAL RESULTS 

We calculate the phase-synchronization index matrix 
of all the sections which consist of continuously recorded 
multi-channel EEG time series. An average is performed 
over all twenty-five sections. Then the brain functional 
networks are generated using an alterable N (the num- 
ber of nearest neighbors). The average cluster coefficient 
(CC) and character path length (CPL) are computed as 
a function of N (see in Fig. Q]). The brain functional 
networks behave large CC and short CPL, which sug- 
gests the small-world property in consistency with previ- 
ous works. 

Community structure and modular architecture are 
two crucial properties of brain functional networks dur- 
ing visual task. The fast GN algorithm is applied to ex- 
plore these hierarchical organizations of these networks. 
Results are shown in Fig. [5] for different N, of which 
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FIG. 3: (Color online) The distributions of (C) and (k„ n ) as 
a function of k. Network is generated with N = 3. 



the minimum value is 3. Figure [2] concretely shows that 
maximum modularities monotonously decrease with TV, 
all above 0.5 even when N < 30, which obviously dif- 
fers from that of randomly connected network. Fur- 
thermore, small maximum modularity also corresponds 
to large mean degree of network, which suggests that 
the more edges bridge the communities to reduce gaps 
of clusters and render community structure less visible. 
Obviously, the number of community corresponding to 
maximum modularity is also a monotonously decreasing 
function of N. With these considerations in mind, we 
mainly investigate the community structure and modu- 
lar architecture of brain functional network with N = 3. 

To better understand the architecture structure of 
brain functional network, we study the average cluster 
coefficient (C) for vertices with degree k. Their relation 
reveals a negative correlation, which suggests that low- 
degree vertices generally belong to well connected clus- 
ters while the neighbors of high-degree vertices belong to 
many different communities which are not directly con- 
nected among themselves, namely hierarchical organiza- 
tion (see Fig. [3]). Simultaneously, the assortative mixing 
pattern (degree-degree correlation of vertices) is inves- 
tigated by using a measure of average nearest-neighbor 
degree (K nn ) that defined as the average over vertices 
with degree k (see Fig. [3]). Degree mixing can be or- 
ganized into two patterns: assortative behavior if (K nn ) 
increases with k, which indicates that high-degree ver- 
tices are preferentially connected with other high-degree 
vertices, and disassortative behavior if (K nn ) decreases 
with k, which denotes that links are more easily built 
between high-degree vertices and small ones. In Fig. [31 
we find that (K nn ) decreases with k, which indicates a 
disassortative behavior of brain functional network. 

The community structure of brain functional network 
is further described by a dendrogram plot computed with 
fast GN algorithms, as shown in Fig0] In Fig. 
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FIG. 4: (Color online) Modularity (a) and dendrogram plot 
(b) when network is generated with N = 3 



ularity is 20, which is marked by a dash-line. Thus, in 
Fig. |4(b)| the dendrogram plot only shows hierarchical 
tree that splits the network into 20 communities. Note 
that we do not present the whole dendrogram plot so that 
the end points of hierarchical tree still denote communi- 
ties with numeric labels from 1 to 20 that are randomly 
ordered in Fig. |4(b)[ 

To further investigate the coincidence of the two kinds 
of hierarchical organizations, we compute the overlap be- 
tween modular architecture of brain functional network 
and functional groups based on anatomy parcellation of 
brain cortex in two ways. On one hand, we check how 
many EEG channels of same community belong to the 
same function group, which is represented by the com- 
position of community written as 



c (i,j) = \c(i)nG(j)\/\c(i)\. 



(3) 



On the other hand, we investigate whether EEG channels 
of the same functional group are exactly divided into an 
identical community, which is described by participation 
of functional group defined as 



P (ij) = \\c(i)nG(j)\/G(j)\. 



(4) 



4(a) the 



number of community corresponding to maximum mod- 



In Eqs. [3] and [H C(i) and G(j) denote the chan- 
nel set of community i and functional group j, respec- 
tively. If there is a perfect one-to-one correspondence be- 
tween communities and functional groups, we will have 
c(i,i) — p(i,i) = 1 and c(i,j) = p(i,j) = 0. Figure [S] 
shows the composition of the communities, which sug- 
gests that most of communities are mainly formed by 
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channel sets that are restricted to the same Brodmann 
area. For instance, the community 12 (numeric label) 
almost overlaps with Brodmann area 9, i.e. the chan- 
nels belonging to the same function group (by anatomi- 
cal parcellation) are densely connected, or coupled in the 
corresponding brain functional network. In Fig. |5(b)[ 
the participation of functional group reveals that most 
of channel sets restricted to the same Brodmann areas 
are divided into the same communities. In particular, 
Brodmann areas 10, 17, and 47 are completely over- 
lapped with the communities 20, 10, and 1 (numeric la- 
bels), respectively. It is noticed that some channel sets 
of functional groups are divided into several large com- 
munities because Brodmann areas involved with these 
functional groups are distributed in both left and right 
cerebral hemisphere with much larger spatial distance 
which reduces the coupling strength among neurons. For 
instance, Brodmann area 21 is mainly divided into two 
larger communities with numeric labels 1 and 15, which 
are distributed in left and right cerebral hemispheres, re- 
spectively. This explains why the one-to-one correspon- 
dence between communities and functional groups are 
not perfect (i.e. c(i,i),p(i,i) ^ 1). 



In a higher level of hierarchy, the coincidence of modu- 
lar architectures (i.e. larger communities) of brain func- 
tional network and functional substructures (see Tab HI)) 
of brain cortex is also studied in aforementioned ways. 
Concretely, the brain functional network is partitioned 
into 9 larger communities. The composition of the com- 
munities shows that the communities are mainly formed 
by channel sets associated with functional substructures 
like Vision, Sensory, Vision-temporal, Motor, which is 
well consistent with the visual task involved in visual, 
judgment and motor functions of brain (see Fig 
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For the participation of functional substructure, most of 
functional substructures are almost divided into unique 
communities, except for the Visual-temporal part, which 
is symmetrically distributed in two hemispheres and is 
separated into two communities, see Fig. 



6(b) 



In addition, the electrical activity of neuronal net- 
works is known to oscillate at various frequencies and 
amplitudes. The densely interconnected oscillators in 
same communities will synchronize more easily than 
those in different communities, which suggests that the 
synchronization reveals hierarchical organization for a 
network with a nontrivial community structure [39TJ43j |. 
Moreover, the simulation of neuronal activity based on 
anatomical cat brain network determines that the cor- 
related clusters are consistent with anatomical areas of 
same brain functions [44Tj46j . Thus, the coincidence of 
modular architectures of brain functional network and 
function groups of brain cortex at two levels of hierar- 
chy suggests that the correlation of EEG time series in 
same functional groups is much stronger than that in 
different ones and the hierarchical organization of brain 
functional network may be a consequence of functional 
segmentation of brain cortex. 
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FIG. 5: (Color online) (a) Composition of community from 
different function groups restricted to Brodmann areas, (b) 
Participation of functional group in communities of brain 
functional network. The result reveals the coincidence be- 
tween modular architectures of brain functional network and 
functional groups based on anatomy parcellation of brain cor- 
tex. Network generated with N = 3 is divided into 20 com- 
munities. 



IV. CONCLUSION 



In conclusion, we have investigated the hierarchical or- 
ganization of brain functional network derived from high- 
resolution synchronous EEG time series during visual 
task through defining vertices as EEG channels and eval- 
uating connectivity between channels using a measure 
of phase synchronization. The resulting brain functional 
networks show common small-world property and com- 
munity structure organized in hierarchical way. Mean- 
while, by clustering EEG channels into functional groups 
based on anatomical parcellation of brain cortex, we find 
that the modular architectures of brain functional net- 
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works are in coincidence with these functional groups at 
different levels of hierarchy via computing overlap be- 
tween communities of channel sets and functional groups. 
These interesting results suggest that population of neu- 
rons performing the same functions excite and inhibit in 
identical rhythms. This is reflected by that the corre- 
lation between pairs of EEG time series (channels) rep- 
resenting the neuroelectrical activity of same functional 
cortex region is enhanced during tasks related to vision, 
judgement and motor functions of brain. The structure- 
function relationship further reveals that the strong con- 
nection of channels and community formation in brain 
functional networks may be a consequence of functional 
segmentation of brain cortex. 
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